(N 

Ph. 



AN EXTENDED CLASS OF ORTHOGONAL 
POLYNOMIALS DEFINED BY A STURM-LIOUVILLE 

PROBLEM 

'o\ 

Q ; DAVID GOMEZ-ULLATE, NIKY KAMRAN, AND ROBERT MILSON 

(N 

P5 ■ Abstract. We present two infinite sequences of polynomial eigen- 

, ^ I functions of a Sturm-Liouville problem. As opposed to the classical 

orthogonal polynomial systems, these sequences start with a poly- 
nomial of degree one. We denote these polynomials as Xi-Jacobi 
and Xi-Laguerre and we prove that they are orthogonal with re- 
spect to a positive definite inner product defined over the the com- 

i^ ■ pact interval [—1, 1] or the half-line [0, oo), respectively, and they 

are a basis of the corresponding L^ Hilbert spaces. Moreover, we 

,S^ • prove a converse statement similar to Bochner's theorem for the 

j^ I classical orthogonal polynomial systems: if a self-adjoint second 

order operator has a complete set of polynomial eigenfunctions 
{Pi}iZij then it must be either the Xi-Jacobi or the ATi-Laguerre 
Sturm-Liouville problem. A Rodrigues-type formula can be de- 

\]^ I rived for both of the Xi polynomial sequences. 

m 

(^ 
en 

O ■ 1. Introduction 

OO ! 

O ■ The classical orthogonal polynomial systems (OPS) of Hermite, La- 

guerre and Jacobi are most often characterized as the polynomial solu- 
tions of a Sturm-Liouville problem, following the celebrated result by 
S. Bochner: if an infinite sequence of polynomials {-Pn(a^)}$^o satisfies 
a second order eigenvalue equation of the form 

p(x)P; + g(x)P^ + r(x)P„(x) = A„P„(x), n = 0,1,2,... (1) 

then p{x),q{x) and r{x) must be polynomials of degree 2, 1 and re- 
spectively [211 E]- In addition, if the {Pn(a^)}^o sequence is an OPS, 
then it has to be (up to an afiine transformation of x) one of the clas- 
sical orthogonal polynomial systems of Jacobi, Laguerre or Hermite 

[D [m El [m [IT] 

Much work has been done since the 1940s until present in different 
generalizations and extensions of these classical families. One main line 
of research has dealt with polynomial sequences defined by differential 
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equations of order higher than two, leading to the Bochner-Krall or- 
thogonal polynomial systems [TB]. For a good review on this subject, 
see for instance |8j. 

When the measure is supported over a discrete set, we speak of dis- 
crete orthogonal polynomials. The equivalent to the classical families 
(Meixner, Hahn, Kravchuk, Charlier, etc.) are orthogonal polynomials 
that satisfy a difference equation of hyper-geometric type instead of a 
differential equation. This topic is reviewed for instance in [T5] . 

Probably the most general class is that of the Askey- Wilson poly- 
nomials, P], a generalization of the classical families that satisfy q- 
difference equations and reduce to the classical and discrete families un- 
der special or limiting cases. Corresponding generalizations of Bochner's 
theorem also exist for polynomials in the Askey- Wilson scheme pT|[T^. 

Another possible generalization concerns the semi- classical orthog- 
onal polynomials, characterized by the fact that the log-derivative of 
the weight factor is a rational function [13]. Semi-classical polynomi- 
als have similar properties as their classical counterparts: they form a 
positive-definite orthogonal family which is complete in the correspond- 
ing L^(ti;)-space, the sequence of their derivatives is not orthogonal but 
quasi-orthogonal [21], they satisfy a second order differential equation 
of the form ([1]), where the coefficients p{x,n), q{x,n) and r{x,n) have 
an explicit dependence on n [H [23]. When the classical weights are 
modified by multiplication by a rational function (with poles and zeros 
outside the interval of orthogonality), the modified weights are semi- 
classical [251 H]- Expressions for these orthogonal polynomials can be 
obtained through the application of Uvarov's determinantal formula 



In the present paper we introduce orthogonal polynomials with ratio- 
nal weights that are eigenfunctions of a Sturm-Liouville operator and 
are therefore fundamentally different from the semi-classical orthogo- 
nal polynomials. The application of Uvarov's determinantal formulas 
gives rise to a sequence of polynomials that begins with the polynomi- 
als of degree zero and consists of polynomials which are orthogonal to 
1. In contrast, the families described below start with a polynomial of 
degree one and are not orthogonal to 1. Our approach leads to novel 
examples that are neither classical nor semi-classical. 

Many of the generalizations referred to above aim to retain the nice 
properties that derive from the Sturm-Liouville character of a classical 
OPS. However, it seems to be a well established fact in the literature 
that no complete orthogonal polynomial systems other than the classi- 
cal ones arise as solutions of a Sturm-Liouville problem. This is indeed 
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the case if the operator belongs to the Bochner class ([I]), as was proved 
by Lesky [IB]. 

We argue that from the point of view of Sturm-Liouville theory this 
restriction is not essential. It has been observed [6] that certain in- 
stances of classical orthogonal polynomial families have the following 
curious property: the polynomials are formal eigenfunctions of the op- 
erator ([1]), but a finite number of initial polynomials are not square 
integrable. Consider, for instance the family of Laguerre polynomials 
P„(x) = L~^{x), n = 0,1,2, .. ., or more specifically 

Po{x) = 1, 

P,{x) = ^x{x-2), 

P2(x) = --x{x^ - 6x + 6), 

D 



P„(x) = xL''^1^{x), n > 1; 

The orthogonality is with respect to the weight W{x)dx = x~^e~^dx, 
which implies that Po{x) is not square integrable. Only the polynomi- 
als Pi,P2,P3, . . . arise as eigenfunctions of the corresponding Sturm- 
Liouville problem, and (therefore) it is this truncated sequence which 
is complete in the L'^(W{x)dx, (0, oo)) space. 
The following question is therefore of interest: 

What sequences of polynomials can arise as eigenfunc- 
tions of a Sturm-Liouville problem? 

The main idea of our paper is to show that the answer to the above 
question takes one outside the realm of classical and semi-classical 
orthogonal polynomials. In other words, if the sequence {Pn}^=m is 
allowed to start with a degree m > 1 polynomial, then there exist 
complete sequences of polynomial eigenfunctions that obey differential 
equations different from ([T]). 

In this paper we treat the case m = 1. In Section [2] we introduce 
the Xi-Laguerre and Xi-Jacobi orthogonal polynomial systems. These 
novel families are crucial to our main result. Theorem 12.11 — a classi- 
fication of complete orthogonal polynomial sequences starting with a 
degree one polynomial that satisfy a Sturm-Liouville problem. This 
theorem can thus be viewed as the corresponding extension of the clas- 
sical results of Bochner and Lesky. 
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Completeness of the new polynomial families is proved in Section O 
using a suitable extension of the Weierstrass approximation theorem. 
Section m completes the proof of the main theorem. Some of the results 
contained in this section rest on the classification of exceptional poly- 
nomial subspaces of co-dimension one and the spaces of second order 
differential operators which leave them invariant. We will use some of 
these results without proof, referring the interested reader to the recent 
publication [TD] where all the details are given. Finally, Sections and El 
describe some properties of the new polynomial families: factorization 
of the second order operator, Rodrigues-type formula, normalization 
constants, relation with the classical families, three-term recurrence 
relation and some basic properties of the zeros. 

By way of conclusion, we mention that since our paper was posted in 
preprint form on the arXiv, the Schrodinger operators and potentials 
for which our new orthogonal polynomials appear as eigenfunctions 
(when multiplied by the corresponding weight) have been determined 
and studied [2U]. The potentials are deformations of the radial os- 
cillator or the Scarf I potential obtained by the addition of rational 
functions, and they are shape invariant. 

2. Definitions and Main results 

Orthogonal polynomial systems are usually understood to start with 
a polynomial of degree 0. However, from the point of view of Sturm- 
Liouville theory, this restriction is unnecessary. The preceding obser- 
vation motivates the following. 

Definition 2.1. We define a polynomial Sturm Liouville problem, or 
PSLP for short, to be a self-adjoint Sturm-Liouville boundary value 
problem with a semi-bounded, pure-point spectrum and polynomial 
eigenfunctions. 

Definition 2.2. For integer fc > 0, we will say that a polynomial se- 
quence {yn}'^=k is degree k {k-V&) if it starts with a polynomial of 
degree k and deg|/„ = n. A fc-PS is a degree k orthogonal polyno- 
mial system (k-OPS) if there exists a positive measure W{x) dx on an 
interval 

/ = (xi, X2), —00 < xi < 0:2 < C)0 (2) 

such that 

(i) The moments are well defined: 

x"iy (x) dx < 00, n = 0,1,2,...; (3) 

I 
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(ii) The polynomials are orthogonal: 

ym{x)yn{x)W{x) dx = 0, m ^ n; (4) 



(iii) The sequence is a basis for the Hilbert space L^(/, W dx). 

We remark that: 

(i) The assumption of self-adjointness in Definition 12. II means that 

the eigenf unctions of a PSLP form an OPS. 
(ii) Item (iii) in Definition 12.21 implies that a 0-OPS cannot be a 

k-OPS for A; > 0. 
(iii) The classical or semi-classical systems are necessarily 0-OPS. 
This means that the fc-OPS represent a genuinely new general- 
ization of the classical OPS, even if they share the property of 
admitting rational weights. 

Definition 2.3. Consider a PSLP whose eigenfunctions form a fc-OPS. 
We call the polynomial system classical if the second-order differential 
equation in question is of Bochner type ([1]). Otherwise we call the 
polynomial system exceptional, or Xk for short. 

It is known that classical Laguerre polynomials with negative integer 
parameters constitute a fc-OPS where k > 1 [6j. Some partial re- 
sults are also available for Jacobi polynomials with negative integer 
parameters [2]. We believe these to be the only classical examples 
where the polynomial eigenfunctions of a second-order Sturm-Liouville 
problem begin with a degree k > 1, but to our best knowledge this 
question has not been explicitly investigated in the literature. Turn- 
ing to exceptional polynomial families, Bochner's result shows that an 
Xq polynomial system is impossible. By contrast, the Xi definition is 
non-vacuous. 

2.1. Xi- Jacobi polynomials. Let a y^ (3 he real parameters and 

«=o(/^-")) b = - , (5a) 

2 p — a 

c = b+l/a. (5b) 
Consider the polynomials 

ui = X — c, Ui = {x — by, i >2, (6) 
the first n of which provide a basis of the space S^''^: 

C' = {peVn\p'{b)+ap{b) = 0} (7) 

= span{Mi,M2,...,Mn} (8) 
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The following restrictions will be required on the real parameters a, /3: 

a > -1, (3 > -1, (9a) 

sgna = sgn/?, (9b) 

the last of which ensures |6| > 1. We define the following measure 

dfla,f3 = Wa,l3dx, XG(-1,1), (10) 

^ _ {1-xni+xf _ {i-xni+xr 

and observe that Wa^p > for —1 < a; < 1 so the scalar product 

{f,9)a,f3-= f{x)g{x)dfia,f3, (12) 



is positive definite. 



>(",/3) 



oo 



Definition 2.4. We define the Xi-Jacobi polynomial sequence < P/' 

L J i=l 

as the polynomials obtained by Gram-Schmidt orthogonalization from 
the sequence {ui}'fZi in & with respect to the scalar product flT^ . and 
by imposing the normalization condition 

From their definition it is obvious that degP„ = n. However, as 
opposed to the ordinary Jacobi polynomials, the sequence starts with 
a degree one polynomial. 

2.2. Xi-Laguerre polynomials. Let k > Ohe a real parameter. Sim- 
ilarly, consider now the sequence 

y^=x + k + l, Vi = {x + k)\ i>2 (14) 

We define the following measure on the interval x G (0, oo): 

dftk = Wkdx, (15) 

and observe that Wk > on the domain in question so the following 
scalar product is positive definite: 

{f,9)k-= f{x)g{x)dfik, (17) 

Jo 
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Definition 2.5. We define the Xi-Laguerre polynomial sequence < L^^ M- 

as the polynomials obtained by Gram-Schmidt orthogonalization from 
the sequence {vi{x)}'^i in (fT^ with respect to the scalar product (TTTj) 
and subject to the normalization condition 

Li^Hx) = -, rr + lower order terms n > 1. (18) 

{n-\)\ 

Note that the Xi-Laguerre polynomial sequence starts with a polyno- 
mial of degree 1. 

Definition 2.6. For a,/? subject to the restrictions ([9]), let 

To^Av) = (^' - 1)/ + 2a (\^) ((^ - ^)y' - y), (19) 

where a, b and c are related to a,/? by ([5]). We define the Xi-Jacobi 
boundary value problem to be the differential equation 

TaAy) = ^y^ (20a) 

where y = y{x) is a twice-differentiable function defined on x G (—1, 1) 
subject to the boundary conditions 

lim (1 - x)'^-^\y{x) - {x - c)y'{x)) = 0, (20b) 

lim (1 + xy+\y{x) - (x - c)y\x)) = 0. (20c) 

X— > — 1 + 

Definition 2.7. For fc > let 

n{y) = -xy" + (^^^ {{x + k + l)y' - y) (21) 

We define the Xi-Laguerre boundary value problem to be the differen- 
tial equation 

Tk{y) = \y, (22a) 

where y = y{x) is a twice differentiable function on x G (0, +oo) subject 
to the boundary conditions 



lim x^+^e"'(y(x) - (x - c)y'(x)) = 0, (22b) 

lim x''^^e'''{y{x) - (x - c)y'{x)) = 0. (22c) 

x— ►oo 

We are now ready to state the main result of this paper. 

Theorem 2.1. The Xi-Jacobi and Xi-Laguerre boundary value prob- 
lems are PSLPs. Their respective eigenfunctions are the Xi-Jacobi and 
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Xi-Laguerre 1-OPSs defined above; we have, 

Ta,pPt^^ = {n-l){a + P + n) P^^), n = 1, 2, . . . , (23) 

r,L(f) = (n-l)LW, n = l,2,... (24) 

Conversely, if all the eigenpolynomials of a PS LP form a 1-OPS, then 
up to an affine transformation of the independent variable, the family 
in question is either a classical 1-OPS, or Xi-Jacobi, or Xi-Laguerre. 

At this point some remarks are due in turn: 

i) Observe that although the components of Ta^ij and Tk are ra- 
tional functions, these operators have an infinite family of poly- 
nomial eigenfunctions. 

ii) Note that both equations belong to the Heine-Stieltjes class 
fr2\ [2B] . i.e. they can be written as py" + qy' + ry = where 
p, q and r are polynomials of degrees 3, 2 and 1 respectively. 

iii) The existence of these new families of polynomial eigenfunctions 
of a second order eigenvalue equation is not in contradiction 
with Bochner's theorem, since one of its premises is that the 
countable sequence of polynomial eigenfunctions should begin 
with a constant. 

iv) Since the sequences start with a first degree polynomial, one 
might think at first that they cannot be dense in the corre- 
sponding L^ space, but we shall see below that this is not the 
case. 

v) The differential expression (IT^ defines a unbounded operator 
on a suitably chosen dense subspace of L^((— 1, 1), rf/iQ,,/?). As 
per the general Sturm-Liouville theory [7], if one takes the max- 
imal such domain and restricts it by imposing boundary condi- 
tions fl20bp (I20cp . one obtain a self-adjoint operator. Alterna- 



tively, one can construct a self-adjoint operator by showing that 
an operator with £^"'* as the domain is essentially self-adjoint. 
This approach is carried out in Section m Similar remarks hold 
for the Laguerre case, 
vi) Both the Xl-Jacobi and the Xl-Laguerre SLPs admit limit- 
point and limit-circle subcases depending on the value of the 
parameters a, f3, k. Details of this analysis can be found in [7]. 

The proof of Theorem 12.11 is based on the classification of Xi sub- 
spaces. For this reason, we recall the necessary results and definitions 
of this classification below, referring the reader to [lOj for further details 
and proofs. 



AN EXTENDED CLASS OF ORTHOGONAL POLYNOMIALS 9 

Definition 2.8. Let M be an n-diinensional subspace of 

Vn = spanjl, X, x^, . . . , x"}. 

We say that M is a codimension one exceptional subspace of Vn {Xi- 
subspace), if there exists a second order differential operator T such 
that TM C M but TP„ ^ P„. 

The main result of the classification of Xi spaces performed in [10] 
states that every Xi-space is projectively equivalent to the space S^''' 
defined in ([7]). For the scope of this study we shall require a stronger 
property, namely that the differential operator T preserves each sub- 
space of the infinite flag 

Sf C ^2'' C ^3'' C ■ • • (25) 

TS^'^ C S^'\ Vn > 1. (26) 

Let a, b, c be real constants related by fl5bl) and set 

pl^x) = k2{x - b)'^ + ki{x - b) + ko, (27a) 

q{x) = a{x-c){ki{x-b) +2ko), (27b) 

r{x) = -a{ki{x-b)+2ko), (27c) 

where ko, ki and k2 are real constants and we assume that k^ 7^ 0. 
Let T define the second-order operator 

T{y):=p{x)y" + ^y' + ^y. (28) 

X — b X — b 

We are now ready to state the following theorem whose proof can be 

found in [TOl: 



Theorem 2.2. The operator T defined in fl28|) with (1271) leaves invari- 
ant S^'^ for all n > 1. Therefore, the eigenvalue equation 

Tyn = KVn (29) 

defines a sequence of polynomials {yn{x)}'i^^i, where yn € £n^ with 
n = degyn and where 

^n = {n — l){nk2 + aki), n > 1. (30) 

Conversely, suppose that T is a second-order differential operator such 
that the eigenvalue equation fl29l) is satisfied by polynomials yn{,x) for 
all degrees n>l, but not for n = 0. Then, up to an additive constant, 
T has the form (1251) subject to flTTI) . and yn € £n^ ■ 

Remark 2.1. The Xi-Jacobi and Xi-Laguerre operators defined in (TT9l) 
and (I2T]) are particular instances of the general Xi operator (128!) with 
(!27|) . In particular, for the Xi-Jacobi take p{x) = x^ — 1 and the 
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parameters a, /5 are related to a, 6, c by (^. For the Xi-Laguerre take 
p{x) = —X and 

a = -l, b=-k, c=-{k+l). (31) 

With the choices above, the general eigenvalue formula (!30!) provides 
the spectrum of the Xi-Jacobi and Xi-Laguerre operators in (!23|) and 



3. Completeness of the Xi-Jacobi and Xi-Laguerre 
polynomial sequences 

In this section we establish the completeness of the Xi-Jacobi and 
Xi-Laguerre polynomial sequences in their corresponding L^ spaces. 
This fact might at first seem counter-intuitive since the classical poly- 
nomial sequences are no longer complete if the constants are removed 
from the sequence. 

Before we prove this result, it is convenient to state the following use- 
ful lemma, essentially a trivial extension of Weierstrass approximation 
theorem, which can also be applied to higher codimension polynomial 
subspaces. 

Lemma 3.1. Let V denote the ring of polynomials in x &M. with real 
coefficients and define V dV to be the following subspace ofV: 

'P= lpeV\J2 <^id'\^i) =0, i = 1, . . . , fc. 
I j=o 

where the k points Xi,...,Xk ^ [~1) 1] and p^^\xi) denotes the j-th 
derivative of p evaluated at Xi. 

Then V is dense in C[—l, 1] with respect to the supremum norm. 

Note that the previous lemma also holds if two or more points Xi are 
allowed to coincide, i.e. if more than one condition is imposed at each 
point. 

Proof. We need to show that given an arbitrary / G C[— 1, 1] and any 
e > 0, there exists a polynomial p &V such that 

\f{x)-p{x)\<e Va;G[-l,l]. 

Consider the function 

g{x) = , ^^""^ G C[-l, 1] 
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since all the poles Xi lie outside the interval [—1, 1]. By the Weierstrass 
approximation theorem, there exists a polynomial p eV such that 

k 

\g{x) - p{x)\ < — Vxg[-1,1], where a = TT(1 + IxiD^^""* 

1=1 

But then, the polynomial p = ni=i(^ ~ Xi)^~^^^ p{x) belongs to V and 
we have 



|/(x) -p{x)\ 



i=l 



iXj Jbrf } 



\g{x) — p{x)\ < e Vxg[— 1,1] 



since \{x - Xi)^+'"»| < (1 + |xi|)^+''' for x G [-1, 1]. D 

Proposition 3.1. // |6| > 1, the space S'^'^ = [j^^n^ ^■^ dense in 

L2([-l,l],#„,^). 

Proof. Since 

S'''^ = {peV\p'{b)+ap{b) =0}, 

and \b\ > 1, Lemma [3.11 ensures that £^°'^ is dense in C[— 1, 1] with re- 
spect to the supremum norm, therefore also dense in L^([— 1, 1], W^^/s)- 

' D 

Proposition 3.2. The Xi-Jacobi polynomial sequence < P>"''^' I is 
a I- OPS. 

Proof. The sequence < P^' > is orthogonal by construction, it suf- 
fices then to prove that it is a basis of L^([— 1, 1], Wa^fs)- But by def- 
inition span < p}-°'''^' I = g"-'^^ so Proposition 13.11 states the desired 
result. D 

In order to prove that the Xi-Laguerre polynomials {L\ }°Zi are an 
orthogonal basis of L'^{[0,oo), fik) we cannot use Lemma 13.1 1 since it 
only applies to the compact case. To this end, we state and prove the 
following: 

Lemma 3.2. The vector space 

£ = {peV\p{-k) = 0} 

is dense on the Hilbert space L^([0, oo),/ifc) where 

d^k = {x + kYdfik = x^e~^ dx. 
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Proof. Since V = M.(BS it suffices to show that 1 is in the L^(/ifc)-closure 
of S. To that end define the function 

,, , fo if 0< X < A;, 

I 1/x it X > k. 

which is clearly in -^^([0, oo), /i^). Since the associated Laguerre polyno- 
mials are dense in L'^{[0, oo), fik), [2Z1, there exists a polynomial p eV 
such that 

POO 

/ |/(x) - p(x) |2a;'=+2e- dx < e'^^e, 
for a given e > 0. Hence, 

/ \l-{x + k) p(a; + /c) px^'e"^ dx = 
Jo 

/■oo 

\l/{x + k)- p{x + k) |^(x + fc) Ve""^ dx 



'0 

POO 

< / \l/{x + k)-p{x + k)\'^{x + k)''+\-''dx 

Jo 

POO 

= e^ \l/x-p{x)\'^x''^'^e~''dx 

Jk 

"OO 



Jo 
< e 

D 

We can now prove the following 

Proposition 3.3. The Xi-Laguerre polynomials \Li > are an or- 
thogonal basis of L"^ {[0,oo), fik). 

Proof. Since < L^ !■ are defined by Gram-Schmidt orthogonalization 

from the sequence {vi{x)}°Zi, the set is orthogonal by construction and 
it suffices then to prove that 

£~^~^ := spa,n{vi}°l^ is dense in L^([0, oo),/ifc). 

Given an arbitrary / G L^([0, oo),/ifc) and e > 0, set 

f{x) = f{x)/{x + k), x>0, 
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and note that / G -^^([0, oo), /i^)- Lemma [3^2] ensures that a polynomial 
p{x) exists such that 

POO 

/ \fix) - {x + k)p{x)\'^x''e~'' dx < e, 
Jo 

Therefore 

/ \fix)-ix + k)Mx)\''dixk<e 
Jo 

and since {x + k^plx) G S^^'^^ this completes the proof. D 

4. Proof of Theorem 12.11 

We begin by recalling some basic facts of Sturm-Liouville theory. An 
arbitrary second order eigenvalue equation 

T{y) = p{x)y" + q{x)y' + r{x)y, 

ny) = Xy, 

can be written in self-adjoint form 

{{pWy')'{x) +r{x)W{x)y{x) = XW{x)y{x), 

provided the function W{x) satisfies a Pearson's type first order equa- 
tion 

{p{x)W{x)y - q{x)W{x) = 0. (32) 

which determines W{x) uniquely up to a multiplicative factor as 

Wix)=pix)-'exp (^I'^de) . (33) 

The following well-known identity establishes the formal self-adjointness 
of T relative to the measure W{x)dx: 

\T{f)g-T{g)f){x)W{x)dx= (34) 



p{x)W{x){f'{x)g{x) - f{x)g'{x)) 



X2 



where — oo < zi < X2 < cc and f{x),g{x) sufficiently differentiable 
functions. The operator T is symmetric if boundary conditions are 
imposed such that the right hand side of (1M|) vanishes. If j/i, 1/2 satisfy 
the eigenvalue equation 

Tyi = Xiyi, i = 1,2 

with Ai 7^ A2 and T is symmetric, we have 

(A1-A2)/ yi{x)y2ix)W{x)dx = 0. 



Xl 
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SO Hi, 1/2 are orthogonal relative to W{x) dx. 

Remark 4.1. The weight function Wa,i3 defined in (1 11 1) satisfies Pear- 
son's equation (!32|) for T = Ta^f^ shown in (fTOl) . Similarly, the weight 
function Wk defined in ( IT6l) satisfies ( l32l) for T = T^ defined in ( 12T1) . 



Forward statement of Theorem 12. IL We can now prove the for- 
ward implication of Theorem 12.11 namely that the Xi- Jacobi and Xi- 
Laguerre SLPs defined in (12(1 and (12^ have a simple, pure-point spec- 
trum bounded from below and a 1-PS of eigenf unctions. 

Let us argue the Xi-Jacobi case and observe that the same argu- 
ments apply mutatis m,utandis to the Xi-Laguerre case. Consider the 
operator Ta^p in (IT9|) defined on the domain £"-'^ . We will show that 
Tq,/3 is essentially self-adjoint. By Theorem 12.21 and Remark 12.11 there 
exist polynomial eigenfunctions ?/„ G £n^ ■, n>l, for which (I20al) holds. 
Since ?/„ satisfy the boundary conditions of the SLP (I2U]) . Remark Wl\ 
and Green's identity (151|) imply that {yn\'^=i are orthogonal with re- 
spect to the weight W^^p in (TTT]) . Moreover, T^^^ is a symmetric and 
semi-bounded operator, so it must have a self-adjoint extension Ta^p 
(see section X.3 in [22]). All the eigenfunctions Hn oiT^^p are also eigen- 
functions of Ta^p and Proposition 13.11 states that {yn}'^=i is a basis of 
L^([— 1, 1], Wa^p). Therefore, the resolution of the id entity associated 
to Ta^/B contains an infinite sum over the corresponding projectors, and 
we conclude that the spectrum is discrete and bounded from below, 
and the self-adjoint extension T^^^ is unique. The spectrum is actually 
given by (12^ . 

In order to prove that the polynomial eigenfunctions {?/„}5^i are 
indeed the Xi-Jacobi polynomials it is enough to note that both se- 
quences span the same fiag of subspaces 

span{yi, ...,yn} = span{pf '^\ . . . , P^"'^)} = £:„"■^ \fn > 1 

and they are orthogonal with respect to the same weight, so up to a 
multiplicative factor they must coincide. 

Converse statement of Theorem 12.11 By assumption T is a second 
order differential operator with a complete set {yn}'^=i of polynomial 
eigenfunctions. Without loss of generality we can assume that the 
coefficients p{x), q{x) and r{x) are rational functions (see Proposition 
3.1 in [TO]). 

If a constant yo is a formal eigenfunction of T, then Bochner's theo- 
rem implies that {yn}'^=o is one of the classical orthogonal polynomial 
systems. Hence, the sequence starting with yi is either not dense (con- 
trary to the assumptions) or constitutes a classical 1-OPS. 
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Let us therefore assume that T has polynomial eigenfunctions for all 
degrees n > 1 but not for n = 0. The converse statement of Theorem 
2.21 asserts that T must be of the form (!28l) with (p7|) . Up to an afiine 
transformation of x, p{x) assumes one of the following five canonical 
forms: 

(35a) 
(35b) 
(35c) 
(35d) 
(35e) 

Writing each of the above 5 cases in self-adjoint form, we obtain the 
following expressions for the weight factor determined by ( l33l) 

(x - ly+'^'^ix + l)"^"'' 



i] 


p{x) = 


= 1-x^, 


u 


p{x) = 


= l+x^ 


iii~ 


p{x) = 




iv^ 


p{x) = 


= X, 


v^ 


p{x) = 


= 1, 



p2a arctan X / 1 i ^2\ab 

u) W{x) = ^^_^y ' (36b) 

m) W{x) = ^^-^ (36c) 

2ax 

Note that the interval cannot include x = b since all eigenpolynomials 
must be square-integrable. We can then use Green's identity flMl) where 
/ and g are any linear combination of the polynomial eigenfuncions of 
T. Theorem 12.21 states that the eigenpolynomials span S"'"'', and since 

{x - hfV C S"^^ 

then Green's identity ([31]) holds, in particular, for 

f[x) = {x- hffi{x), g{x) = {x- hfgi{x), 

where fi.,gi are arbitrary polynomials. We observe that 

f'{x)g{x) - f{x)g\x) = (x - hY{f[{x)gi{x) - fi{x)g[{x)) 

= (x — byh{x), 

where 

h{x) = f[{x)gi{x) - fi{x)g[{x), 
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is an arbitrary polynomial. Since the left-hand side of (j3^ vanishes by 
assumption, the expression 

p{x)W{x){x-b)^h{x) =0 (37) 

J XI 

must vanish for all polynomials h{x), which implies that 

{pW){x^) = ipW){x2) = 0, (38) 

where the above evaluations have to be understood in the limit sense 
if one or both of the endpoints Xi, X2 are infinite. 

It is clear that condition (!38l) excludes cases (ii) and (v) for all possi- 
ble choices of xi, X2- Case (iii) is also excluded by the requirement that 
all eigenpolynomials be square-integrable relative to Wdx on [xi,X2]. 

Case (i) leads naturally to the Xi-Jacobi SLP. Equation (!38l) implies 

Xi = — 1, X2 = 1, ab±a>—l. (39) 

Setting 

a = ab — a, P = ab + a, 

we obtain fl5aj) and the conditions on a, (3 given at the beginning of 
Section [2l In particular, equation ( !39l) implies ( !9al) while ( l9bl) has to 
be imposed to ensure that b lies outside the interval [—1, 1]. With these 
restrictions, the weight (I36al) specializes to the Xi-Jacobi weight Wa,i3 
shown in fITT]) . Theorem 12.21 implies that the eigenpolynomials of the 
SLP are the Xi-Jacobi polynomials. 

Similarly, case (iv) corresponds to the Xi-Laguerre SLP . By rescal- 
ing X we can assume that a = — 1 without loss of generality. The 
condition (155]) implies then 

xi = 0, a;2 = +00, 6 < 1. (40) 

However, for b to lie outside [xi,X2], we must impose 6 < 0. Setting 

a = —1, b = —k 



-) 



we obtain the Xi-Laguerre weight (TT6!) by specializing the weight shown 
in (]36d|) . The same argument as above shows that the given SLP has 



to be the Xi-Laguerre SLP. D 
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5. Properties of Xi-Jacobi polynomials 

5.1. Factorization and Rodrigues formula. Define the following 
lowering and raising operators: 

= (^H^Af^y (42) 

X — b dx \x — cj 

BaAy) = i^^ - 1) (^^] {y + ay) - a{x^ - 2bx + l)y (43) 

\x — cJ 

where a^h^c are related to a, /5 by ([5]), and where the weight W^q,/? is 
defined in (fTTI) . Using the above operators we can factorize Tq^/^ in two 
different ways: 

Tafi = Ba,pAa,p (45a) 

= Aa-1,13-1 Ba-i,p-i -a- p. (45b) 

Another consequence of (HTj) (H51) is the following adjoint-type relation 

{Aa,f3 /, g)a+l,f3+l (/, 5a,/3 5')a,/3- (46) 

relative to the the inner products defined in (TTT]) (TT2!) . 

By virtue of the intertwining relations (HSl) . the above raising opera- 
tor can be applied iteratively to construct the Xi-Jacobi polynomials. 
The difference with respect to the classical raising operators is that on 
each iteration a parameter needs to be shifted by an additive constant. 
More specifically, the following relations hold: 

Aa,pPt''^l{n + a + 13) PtV^'^'^ (47) 



B^,,Pt^'^^^'^ = 2nP!f^P (4J 

I 
For fixed a, P, set 



where Pn is the n-th Xi-Jacobi polynomial. 



6, = 


= b + 


]/a, 


W, - 


= Wa+j,l3+j, 


(y) - 


d 
dx 


r (x-6,)^ 


_{x - bj-i){x - bj+i) 



(49) 
(50) 

(51) 
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Iterating fHSl) and using the identity 

we obtain the following Rodrigues-type formula for the Xi-Jacobi poly- 
nomials: 

(-2)»(„ - l)!Pr' = W--^--)[("-''")^^'^'-] (52) 

(x - hi)WQ 

5.2. Norms. The square of the norm of the Xi-Jacobi polynomials is 
given by 

where 

"~ a + /5 + 2n + l r(n + l)r(a + /3 + n + l)' ^ ' 

The above should be contrasted with the norm formula for the classical 
Jacobi polynomials, namely: 

(1 _ a;)"(i + a;)/3 (P^"'/^))' dx = Cn. (55) 

-1 

5.3. Relation to classical polynomials. The Xi-Jacobi polynomi- 
als Pn' are related to their classical counterparts Pn' by the fol- 
lowing 3-term linear combination: 

pi-^P) = _i(a; _ 5)P^"f + iili^ i^, (56) 

2^ ' ""-^ (a + /3 + 2n-2)' ^ ' 

where b is given by ([S]). Using the 3-term recurrence relation for the 
classical Jacobi polynomials, relation ( l56l) can be rewritten as 

Pi"'^) = -fn P^C''^ + 26^. Pi^f ) - /.„ Pi!'? , (57) 

where 

r ^ n{a + P + n) 

^" (a + /3 + 2n-l)(a + /5 + 2n)' ^ ^ 

= (a + n)(/3 + n) 

^" (a + /5 + 2n-2)(a + /5 + 2n)' ^ ^ 

, ^ (a + n)(/3 + n) 

(« + /? + 2n-2)(« + /3 + 2n-l)' ^ ' 

and where 6 is given by ([5]). Relation (1571) can be established by means 
of (147|) and by using the series definition of the classical Jacobi poly- 
nomials. The details are left to the reader. Using (l53l) (155|) and the 
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orthogonality properties of Pn \Pn' , relation (1571) can be inverted 
to obtain the following identity: 

-\{^- hf pi"''^) = /n+iPi:f ^ - 2 6^„pi:f + kpic''\ (61) 



where 



^ ^ {n + a){n + (5) 

^" (a + /5 + 2n)(« + /? + 2n + 2)' ^ ' 

h = {n-l + a){n-l + P) 

{a + 13 + 2n){a + 13 + 2n + l)' ^ ' 



5.4. Recursion formula. Using (l57j) and the classical 3-terni recur- 
rence identity we obtain the following expression for the classical Jacobi 
polynomial in terms of its Xi counterparts: 



(&'-l)Pi"'^^ = (« + ^)(/3 + ^)(-/n-MPi?2 + |^S) (64) 

- 2(a2-l)6^„Pi:f 

- {a + n + l){l3 + n + l)KPi"'\ 



where a, 6 are given by ([5]). Combining the above identity with fl6Tl) 
yields the following 3-term recurrence for the Xi- Jacobi polynomials: 

/„+i[(62-l)-(a + n)(/5 + n)(a;-6)2]Pi:f + (65) 

-2hg^[{h' - 1) + (a^ - l)(x - hf]Pi^;^'^ + 

^(« + n)(/? + n)x(x-6)2pi^f + 

hn[(h''-l)-{a + n + l){P + n+l){x-hf]P^^'^^ = 0. 

5.5. First few Xi-Jacobi polynomials. The first few Pn polyno- 
mials are: 

p(a,/3) ^ _1 2+a+P 

-^1 2"^ 2(a-/3)' 

p(a,/3) _ a+/3+2 2 _ a^+/3^+2(a+/3) _ a+/3+2 
^2 4 -^ 2(a-/3) -^ 4 ' 

p(a,/3) _ (a+/3+3)(a+/3+4) 3 _ (3+a+/3)(6a+3«^+6^-2a/3+3/3^) 3 

(9«+3a^+9/3+2a/3+3/32) _6a+a2+a3_6^_ga/j_a2/3+/32_a/32+/33 

16 ^ 16(a-/3) ' 
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5.6. Zeroes of Xi-Jacobi polynomials. Many properties of the ze- 
ros of the Xi-Jacobi polynomials follow from the fact that they are 
eigenfunctions of a Sturm-Liouville problem. However, we choose to 
give a direct proof below independent of Sturm-Liouville theory. 

Proposition 5.1. Assume without loss of generality that a < 0, then 
the n-th Jacobi polynomial Pn' (x) has one zero in (— oo, b) and n — 1 
zeroes in f— 1, 1). 



Before proving Proposition 15. 11 let us state the following two lemmas: 

Lemma 5.1. Let P G S^'^ be a polynomial with n real roots. If a < 
and P{b) ^ 0, at least one of these roots lies in (— cxd, b). 

Proof. From ([7]) it follows that 

a = -P'{b)/P{b), 

hence P{b) and P'{b) have the same sign. By Sturm's root counting 
theorem, It is clear that a root of P has to lie in (— oo, b) otherwise P 
cannot have n real roots. D 

Lemma 5.2. P!C''^\b) y^ 0. 

Proof. First note that the Pn' (x) are defined recursively by (l48ll and 
( 144|) . Using ([n]) in (jHD it is clear that (glD has the form 



B.,,y=i^-b)f{x)-^-^—^^y 

where f{b) ^ and g{b) ^ 0. Since P|"' (6) ^ 0, it follows by induc- 
tion that Pt''^\b) ^ for all n > 1. 

D 

Proof of Proposition \5.1[ From the two previous lemmas it follows that 
Pn ' (x) has at most n — 1 zeroes in (6, oo), and in particular at most 
n — 1 zeroes in (—1, 1). Suppose that Pn (x) has ^i, . . . , ^fc, 1 < fc < 
n — 2 zeroes in (—1, 1), and let 

Qi{x) := {x-^i)---{x-^k)- 

If Pn (x) has no zeros in (—1, 1) then take Qi{x) = 1. By Lemma 



5.11 the polynomial Qi ^ £^"'^(x) but we can always choose ^ so that 

Q{x) := {x - OQiix) e S:,^\. 
This is clear because imposing ([7]) on the above expression leads to 
{b-OiQ'iib) + aQiib)) + Qiib) = 0, 
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which can be solved for ^. Again, Lemma [5.11 imphes that ^ ^ (— 1, 1), 
and therefore the function Q{x)Pn (x) does not change sign for x G 
[—1, 1]. Hence 

\ / a,l3 

but this is impossible since Pn is orthogonal to ^^Li- We conclude 

then that Pn has exactly n—l roots in (—1, 1). The remaining root 
has to be real and Lemma [5.11 implies that it lies in (— oo, 6). D 



6. Properties of Xi-Laguerre polynomials 

6.1. Factorization and Rodrigues formula. Define the following 
lowering and raising operators: 



My) 



jx + k + l) 

(x + k) 

{x + k + lf d 



[y' -y) -y, 
y 



Bk{y) 



X ■ 



x + k 
(x + k) 
'x + k + ] 



dx 



x + k + 1 



iv' -y) + ky 



{{x + k + l)m)-'^ 



{x + k + iy 

x + k 



Wk+iy 



(67) 
(68) 
(69) 
(70) 



where the weight W^ is defined in (1161) . Note that we can factorize the 
second-order operator in fl2T|) in two different ways: 



n 



BkAk 
^fc-i-Bfe-i 



(71a) 
(71b) 



and observe the following relations relative to the the inner products 
defined in (|T6l) (|T7|) : 



{Akf,g)k+i = {f,Bkg)k- 



(72) 

By virtue of the intertwining relations flTTl) . the raising operator Bk 
can be applied iteratively to construct the Xi-Laguerre polynomials. 
The difference with respect to the classical raising operators is that on 
each iteration a parameter needs to be shifted by an additive constant. 
More specifically, the following relations hold: 

■-{k+i) 



AkW = L 



n-l 



BkW^^^ =nL, 



{k) 

n+l 



(73) 
(74) 
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where Ln is the n-ih. Xi-Laguerre polynomial. Iterating (1741) we obtain 
(n-l)!LW = (5fe---5,+„_2)^f+"-^\ n = 2,3,... (75) 
Fix k and set 

Using fl7U]l and 



(x + fc + j)^ 

_{x + k + ] -l){x + k + i + l)^ 



(76) 



Lf) = -(a; + A; + l), 

we rewrite fl75l) to obtain the following Rodrigues-type formula for the 
Xi-Laguerre polynomials: 

{x + k + l)Wk 

6.2. Norms. The square of the norm of the Xi-Laguerre polynomials 
is given by 

(i<'))^d.= <^^±^/C.. (78) 

The above relation follows by induction from (1211) 071ap (172]) (175]) . Con- 
trast the above to the norm formula for the classical Laguerre polyno- 
mials, namely: 



oo 



3:V-(4^ydx ^^^ + f + ^^ ^ ir„ (79) 

n\ 

6.3. Relation to classical polynomials. The Xi-Laguerre polyno- 
mials Ln are related to the classical Laguerre polynomials Ln by the 
following simple relation: 

^i'^ = -(^ + A: + l)4'^i + 4'i- (80) 

Using the 3-term recurrence relation for the classical Ln , 

nL(f) + {x-2n~k + l)4''ii + {n + k- 1)4-2 = 0, (81) 
relation (I5U]) may be rewritten as 

4'^ = ^4'^ - 2(n + fc)4'=], + (n + fc)4'=]2- (82) 

The above identity follows by induction from ( 1731) and from the follow- 
ing properties of the classical Laguerre polynomials: 

r(fc) _ r(fc+l) _ r(fc+l) (Qn\ 
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Using (1781) (1791) and the orthogonality properties of Ln and Ln we 
can invert ( l82l) to obtain the following identity: 

(x + k^L^:^ = in + l)Li% - 2{n + A;)lSi + {n + k - l)!^. (85) 

6.4. Recursion formula. Identities ( 1811) and (15^ imply the following: 

(^ + l)(^ + A;)Li% + (86) 

+ (n + A;)(x-2n-A;-l)Ljf_Ji + 

+ (n + A;-l)(n + A; + l)L(f) = fcL^f). 

Combining this with (l85l) yields the following 3-term recurrence for the 
Xi-Laguerre polynomials: 

{n + l)[{x + kf{n + k)-k]L^;:l^ + (87) 

+ {n + fc)[(s + kf{z -2n-k-l) + 2k] L^^l^ + 

+ {n + k - l)[{x + k)\n + k + 1) - k] lI^^ = 0. 

6.5. Zeroes of Xi-Laguerre polynomials. 

Proposition 6.1. The n-th Laguerre polynomial Ln (x) has one zero 
in (— C)0, —A;) and n — 1 zeroes in [0, oo). 

Proof. The proof follows the same arguments as Proposition 16.11 D 

" (k) 

6.6. First few Xi-Laguerre polynomials. The first few Ln poly- 
nomials are: 

Lf^ = -x-{l + k), (88a) 

Lf = x^-k{k + 2), (88b) 

Lf = -lx' + k±s^^ + '<t+^x-l{3 + Ak + e). (88c) 
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